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Abstract 

The interest in the properties of quantum systems, whose classical dynamics are chaotic, 
derives from their abundance in nature. The spectrum of such systems can be related, in the 
semiclassical approximation, to the unstable classical periodic orbits, through Gutzwiller's 
trace formula. The class of systems studied in this work, tiling billiards on the pseudo-sphere, 
is special in this correspondence being exact, via Selberg's trace formula. 

In this work, an exact expression for Green's function and the eigenfunctions of tiling 
billiards on the pseudo-sphere, whose classical dynamics are chaotic, is derived. Green's 
function is shown to be equal to the quotient of two infinite sums over periodic orbits, where 
the denominator is the spectral determinant. Such a result is known to be true for typical 
chaotic systems, in the leading semiclassical approximation. From the exact expression for 
Green's function, individual eigenfunctions can be identified. 

In order to obtain a semiclassical approximation by finite series for the infinite sums 
encountered, resummation by analytic continuation in h was performed. The result, a semi- 
classical approximation, is similar to known results for eigenfunctions of typical chaotic 
systems. The lowest eigenfunctions of the Hamiltonian were calculated with the help of the 
resulting formulae, and compared with exact numerical results. A search for scars with the 
help of analytical and numerical methods failed to find evidence for their existence. 
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1 Introduction 



The spectrum of chaotic systems was studied extensively in recent years. Most of these studies are 
in the framework of the Gutzwiller trace formula ^ for the density of states, that holds in the 
semiclassical limit, and was used successfully to analyze spectral correlations [|3[0,|5|]. The main 
problem with it's application to the calculation of energy levels, is the fact that the Gutzwiller trace 
formula involves a non-absolutely convergent series. Resummation techniques enable to approximate 
the infinite sum by a finite one, and to obtain the spectrum in the semiclassical approximation 
For billiards on surfaces of constant negative curvature, the density of states is given by the exact 
Selberg trace formula, that corresponds to the Gutzwiller trace formula. Since it is exact, it enables to 
distinguish between consequences of the semiclassical approximation of the Gutzwiller trace formula 
and those of other approximations, that are usually introduced, like the diagonal approximation or 
resummation. 

Exploration of the eigenfunctions followed the studies of the spectrum. It is important because their 
properties determine matrix elements which control the nature of important physical phenomena, such 
as transport, dissociation of molecules, absorption and other physical processes where cross sections 
are important. Since the classical dynamics of many real systems are chaotic, the eigenfunctions of 
such systems are of special interest. Besides the physical motivation to study such eigenfunctions, 
they are also very interesting as they exhibit some surprising features. Prominent among these are 
scars, discovered by Heller |^] in numerical calculations. A scarred eigenfunction describes a state of 
the Hamiltonian of the system, whose probability density is unusually large in a narrow tube around 
an unstable classical periodic orbit. This effect was one of the factors that spurred the study of 
eigenfunctions of chaotic systems, which aimed at finding the explicit role classical periodic orbits 
play. Another important factor are mathematical theorems that state that the measure of scarred 
eigenfunctions (out of all eigenfunctions) is zero. 

In the present paper, an exact formula of eigenfunctions for some billiards on surfaces of constant 
negative curvature will be derived. It's main importance is for the understanding of the various 
approximations that were made in order to obtain the corresponding results in the semiclassical 
approximation for more realistic model systems. 



Bogomolny studied the semiclassical limit of Green's function of chaotic systems. He showed that 
it has two kinds of contributions, the first being that of "zero length orbits" and the second consisting 
of a sum over all the classical periodic orbits. From this expression he obtained an expression for 
smoothed eigenfunctions, in which the smoothing is done over a small interval of energy. Because 
of this smoothing the study of individual eigenfunctions was prevented. Later Berry [|ll|] revealed 
the role periodic orbits play in the semiclassical limit of the spectral Wigner function. This function 



is closely related to Green's function, and is very suitable for semiclassical studies [12|. Individual 



eigenfunctions are related to the poles of the spectral Wigner function. Agam and Fishman [13| cast 
the spectral Wigner function in a form that reveals it's poles. They showed that it can be written. 
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in the semiclassical limit, as a ratio between two infinite sums over the classical periodic orbits. The 
denominator is the spectral determinant, whose zeros are the eigenenergies. This enables easy singling 
out of a specific eigenfunction. As the sums that appear in the numerator and the denominator were 
not absolutely convergent, they were treated by a resummation method previously invented by Berry 
and Keating ||7| for the spectral determinant. Later on, Fishman et al. [14| showed that as Green's 



function is related to the solution of a Fredholm integral equation, it is given exactly by a ratio 
between two entire functions. The parameter in which the functions are entire is a book-keeping 
parameter that counts the number of times a surface of section through the system is crossed by 
Feynman paths. In addition, Fishman et al. showed that invoking the semiclassical approximation, 
whose leading order arises from points near periodic orbits, Agam and Fishman's results are obtained. 



The term "scars" was quantified by Agam and Fishman [15| who proposed to study the weight of 
the eigenfunction in a narrow tube around the classical periodic orbit. They found, for some model 
system, that the weight of scars in the eigenfunctions tends to zero in the semiclassical limit of 
vanishing wavelength. This decrease is very slow, in agreement with the observation of scars in 
numerical calculations. 

In the present work, the diagonal part of Green's function of a billiard on a surface of constant 
negative curvature will be expressed as the ratio: 

GMC,C) = ^^, (1.1) 

where A(E') is the spectral determinant (the precise meaning of Ge{CjO w'" elucidated in 
Section Exact expressions for M{CjE) ^^d A(£'), in terms of classical orbits, will be obtained. 
These are similar, in their general structure, to expressions that were obtained in the framework of 
the semiclassical approximation, supplemented by other plausible approximations. It provides further 
evidence for the validity of the approximate formulas that were obtained for realistic systems. 

The special chaotic systems that are investigated are billiards that are the primitive cells of lattices on 
a surface of constant negative curvature. This surface is called the pseudo-sphere. The reason for this 
choice is that the Green function of these systems is given exactly in terms of known contributions 
from classical orbits. In spite of the fact that these systems are chaotic, an exact expression for 
the Green function can be obtained by the method of images, since the surface is perfectly tiled by 
images of the billiards. 

The exploration of scars is of special interest for billiards on surfaces of constant negative curvature. 
Rudnick and Sarnak ||l^ defined scars as eigenfunctions with weight that is divergent in the semi- 
classical limit. They proved that for a sub-class of tiling billiards on the pseudo-sphere ("Arithmetic 
Billiards") there are no scars according to their definition. This is in agreement with the extrap- 
olation of the semiclassical formula of Agam and Fishman. It does not make a clear statement 
about eigenfunctions for small, but finite, wavelengths. This issue was studied numerically for tiling 



billiards on surfaces of constant negative curvature ||17|,|18|,|19[1, and no scars were found. The quest 



for understanding of the difference between systems where scars are found and those where they are 
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not was one of the motivations for developing analytical formulas for eigenfunctions for billiards on 
surfaces of constant negative curvature. 



The outline of the paper follows. In the next two sections, results of previous research are reviewed - 
Section § is a brief summary of previous results concerning the semiclassical behaviour of typical 
chaotic systems. Section ^ introduces the pseudo-sphere and tiling billiards on it, which are the 
specific chaotic systems analyzed in this work. In section Section § the expression (1.1) for Green's 
function is obtained. The exact eigenfunctions densities are it's residues. In order to truncate the 
infinite series that arise. Berry and Keating's resummation method [0] is used in section Section |5|to 
yield a semiclassical approximation to the exact expression. In section Section ^some of the resummed 
expressions are compared to the exact eigenfunctions found numerically. In section Section 0, the 
results of this work are compared with the ones found for general chaotic systems in the semiclassical 
approximation. The existence of scars is studied in section Section B 



Many of the results in the paper rely on the review of Balazs and Voros |^ and on the paper of 
Agam and Fishman [13|. Equations in these papers will be referred to as "BV" or "AF" followed by 



the equation number in the corresponding paper. 



2 Semiclassics 



The starting point of the discussion is the Gutzwiller trace formula, which relates to the density of 
states of typical chaotic systems to the periodic orbits of the corresponding classical system. This 
formula is not convenient, as the density of states is a very singular function of energy. Therefore 
the spectral determinant, whose zeros lie on the spectrum of the Hamiltonian, is presented. In the 
semiclassical approximation the spectral determinant can also be written as an infinite sum over 
periodic orbits, that is obtained from the Gutzwiller series for the density of states. As similar sums 
appear many times in this work. Berry and Keating's method of resumming it to obtain truncated 
series is outlined. Finally, Agam and Fishman's results for eigenfunctions of typical chaotic systems 
are presented. 



2.1 Gutzwiller's trace formula 



The resolvent operator is given by: 

HiE + ie) = ^ ^ , 

E + ie-n 
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where IrL is the Hamiltonian operator. The retarded Green function is the representation of "^(-E+ie) 
in the coordinate-space basis: 



G+(r',r) 



(r'l^lr) = Y: 



C(r')V'o 



E + ie- Ea 



(2.2) 



and the density of states is related to it's trace: 

p{E)=J26{E-E^) = --Q{ 

^ IT ^ 



Tr 



n{E + ie) } 



(2.3) 



-.0+ 



One would like to calculate the spectrum. Very often it is only possible in the leading order of the 
semiclassical approximation. For integrable systems, it can be done by the EBK method |21], that 
relies on the existence of d "good" constants of motion, where d is the number of degrees of freedom. 

Gutzwiller ^ studied p{E) in the semiclassical limit for systems whose classical dynamics are very 
chaotic (hard chaos). In the semiclassical limit, the density of states can be divided into two parts 
with different analytic properties. The first of these is a smooth function of energy and the second 
is an oscillating function of it: 



p{E) ~ p{E) + ppo{E) ash^O . 



(2.4) 



The smooth part originates in the the system's short (non-periodic) classical orbits, and it's leading 
order is given by: 

'^'^5iE-nU,p) 



as h . 



(2.5) 



The oscillating part arises from the classical periodic orbits, and it's leading order is given by the 
Gutzwiller trace formula S: 
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(S'p-j7p)r 



r=l 



as h —>■ 



Det (m: 



(2.6) 



where J2p denotes summation over all the primitive periodic orbits of the system, Tp denotes a 
primitive periodic orbit period, Sp it's action, 7^ the Maslov index associated with it and Mp 
denotes the monodromy matrix for motion perpendicular to it in it's vicinity in phase space. 



As (|2.6| ) stands it does not converge absolutely, due to the exponential proliferation of primitive 
periodic orbits with their length, typical of chaotic systems. Therefore, it is necessary to to treat 
Gutzwiller's trace formula further in order to obtain useful expressions. Such a treatment is presented 
in the next subsection. 
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2.2 Spectral determinant and it's resummation 



For the calculation of the spectrum it is more convenient to study the spectral determinant. The 
semiclassical limit of this function was studied by Berry and Keating [§, ^. in this section some of 
their results will be presented. 



The regularized spectral determinant is defined by: 

A{E)=l[A{E,Ea)iE-Ea) 



(2.7) 



and it is convenient to study because it has zeros on the spectrum of the quantum Hamiltonian, 
instead of the poles the density of states ( ^ ) had. A{E,Ea) is just a regularizing function [22|. 

Berry and Keating showed that the spectral determinant is given asymptotically by a dynamical Zeta 
function [0] (specializing for simplicity to two dimensions): 



A{E) ~ B{E)e 



-i-KN{E) 



nn{i 

p m=0 



,j-Sp-nfo p-(5+m)«s 



I as ?i ^ , 



(21 



where N{E) is the smooth part of the counting function of eigenenergies, given asymptotically Weyl's 



expansion in powers of ?i pl 



are the eigenvalues of Mp, and the uJs are called instability 



exponents. The double product in Eq. can be written as an infinite sum: 

A{E) ~ B{E)e-'^^^'^^ J2 C^,e^^^ as 7i ^ , 



where ji denotes the list of composite periodic orbits ordered by length. In other words: 

= jp^p jp = 0, 1, 2, . . . 



and the coefficients Cn are given by: 



n(-i 



n Det (mi 



(2.10a) 



(2.10b) 



Eq. 2.S suffers from the same convergence problems Gutzwiller's trace formula does. In order to 
treat them Berry and Keating introduced the resummation of (|2.9| ) by analytic continuation from 
the region of complex ^-plane where it converges. The method found to be effective made use of a 
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Cauchy integral, and the result of the analytic continuation was found to be |^: 



A{E) ~ 2^< 



as h ^ 



(2.11) 



where: 



Q\K,h,E) = K' + . 



(2.12a) 
(2.12b) 



This choice ensures 



iiT is a free tuning parameter chosen from the range < K « jpi, 2' 
that number of composite periodic orbits with actions in the cut-off region of the complementary 
error function is much smaller than the number of composite periodic orbits with action less than 
the center of the cut-off as h ^ 0. 



The subsequent terms in the asymptotic expansion for A{E), whose leading term is given by (2.11), 
are given by high order derivatives of the complementary error function. These derivatives are 
proportional to Gaussians centered around the complementary error function cut-off. 



One observes that the resummation causes the effective truncation of the problematic sum in (|2.9|) . 
The center of the cut-off is the composite periodic orbit whose period is half the Heisenberg time M: 



T:{E) = nhp{E) 



(2.13) 



2.3 Eigenfunctions of chaotic systems 

Spectral Wigner's function is a function whose poles lie on the spectrum of the Hamiltonian, and 
whose residues are closely related to it's eigenstates. It's semiclassical behaviour was studied by 
Berry |11,^. Then, Agam and Fishman [13| showed, that in the semiclassical limit, it can be 



written as a ratio between two functions given by sums over periodic orbits. The denominator of 
the expression they found is the spectral determinant of the previous section. This allowed them to 
single out expressions for eigenfunctions. 



Fishman et al. ||1J] were able to put the results of Agam and Fishman |[lj] on a more rigorous footing. 
They showed that the exact Green function is given by a ratio of two entire functions. Then they 
managed to show that in the semiclassical limit their result is identical with Agam and Fishman's 
result. 

The method Fishman et al. employed was Fredholm's method for the solution of certain integral 
equations. In this method the exact solution of the integral equation is given by the ratio between two 
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functions. Each of the functions is given by an absolutely convergent series in powers of a transfer 
operator (T) in a surface of section of the system. Using a semiclassical expression for T, the series 
then assumes the form of absolutely convergent sums over the classical orbits of the system. Then, 
using the same approximations used in [O], results that concurred completely were obtained. 



The expressions found by Agam and Fishman suffer from the same convergence problems Eq. |2_ 
and (|2.9| ) do. This problem was solved for denominator by analytic continuation in the ^-plane by 
Berry and Keating |^], as was mentioned in Section 2^. The same treatment was implemented by 
Agam and Fishman for the numerator. Once this was done, they could present useful expressions 
for eigenfunction densities. The structure of these expressions is reminiscent of Gutzwiller's trace 
formula, for they have a smooth contribution and a periodic orbit contribution. For later reference 
only the latter contribution will be conveyed here (AF 3.64): 



9|| 



as h ^ 



(2.14) 



where q|| is the velocity along the primitive periodic orbit p, X 
the phase space coordinates perpendicular to it, and: 



iQi,---Qd-i,Pi,--- Pd-i) are 



m=0 tJ. 

xe-'"^(^-)+^^^--Erfc 



— XR^X 

h ^ 



CAt,p(^; Eg) + XRpX 



(2.15) 



Confining the discussion to two dimensions, the argument of the complementary error function 
includes the function: 

C,A^,E)=S^,^iE) - ^ , (2.16) 



2'n:h 



and Q'^{K,h, Ea) that was defined in Eq. 2.12b. The coefficients cjf'"^^ arise from expressing (AF 
3.26): 



oo 

Jl (^1 _ etV-*Ve~(Hj)«p'^ X e^^*'-*Tpe"(5+"^)"p , (2.17) 



as a sum over the composite orbits introduced in Section 2.2 



(2.18) 
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where, 

5^,p = Sf, + Sp = {jp + l)5p + . (2.19a) 

S^-p denotes the composite action 5^ after the contribution of the primitive periodic orbit p has 
been removed. The expression for the cjf'"''^ given by (AF C.ll) is identical with: 



j=0 



where C^-p is the coefficient of the composite orbit S^-p (see Eq. |2.19a| ) of (^^), and dj{x) is 
given by (AF C.2): 



dj{x) 



n (1 - 

r=l 



The matrix Rp(x) is related to the eigenvectors of the monodromy matrix: 



Rp(x) 



Rp Rp \ Rp ~^ Rp 



-2R+R^ R+ + R^ 



in which R^ = . 



The function gm [x] appears through an expansion in powers of e (AF 3.11 &c 3.13): 

exp 



/det(M; + /) 
where r was used for repetitions. 



m=0 



h ^ 



In Appendix A it is shown that: 



giXRpX tanh ^ 



xe 



C/j,p(^i -E'en) + XRpX 



(2.19c) 



as /i ^ 



in which another definition was introduced: 

C^P'^') = C^_pe-*(^''^+^)^'^(-l)^'^-^e-5(j'*'-^>^dj-^_j- (e""'') 



(2.20) 



(2.21) 



(2.22) 



(2.23) 
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Finally, the integral in (2.14) can be performed using the stationary phase approximation. The result 
is (AF 4.3): 



E 

p m 



1 IhTT 1 

2i2 



M j=0 



sinh {j + l)n«; 



as 7i ^ . 



(2.24) 



3 Tiling billiards on the pseudo-sphere 



The systems considered in this work are tiling billiards on the pseudo-sphere. There is extensive 
mathematical literature concerning these systems. Various aspects of geometry can be found in |24], 



while many of the results about quantum mechanical properties of these systems, to be conveyed in 
the next subsections, can be found in |||,|2|,|2^. Useful reviews are p0| , |28| 



The pseudo-sphere is a surface of constant negative Gaussian curvature of Setting R = \, all 

distances are measured in units of R. Because the surface is two dimensional, it can be projected 
onto the the complex plane, so that the two dimensional coordinates are given by complex numbers. 
Two very useful projections are the Poincare disk, in which the surface is projected onto the interior 
of the unit disk, and the Poincare half plane, in which it is projected onto the upper half of the 
complex plane. In both of this models the geodesies on the surface turn out to be arcs of circles 
perpendicular to the boundary of the model. 

The invariant measures in both models are: 

, , , 4r drdS ^ , ,■ , -■ x 

dijL{z) = 2 'O'' (3.1a) 

(1 — r^) 

dxdy 

dfJ-iC) = — ^ foi" the half-plane, (3.1b) 

yZ 

where z = re^'f' are disk coordinates, and C, = x+iy are half-plane coordinates. The distance between 
two points is: 

Z \ z — z\ 

cosh 4,2' = 1 + — tTjTT; TTm fo'' the disk, (3.2a) 

(1 — |zp) (1 — |r p) 

1 ic - cP 

cosh dec = 1 + --^^ ^ for the half-plane. (3.2b) 
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There are several kinds of isometries that are going to be encountered in this paper: 



1. Even boosts : these isometries are the analogues of Euclidean translations. They have two fixed 
points at infinity, which lie on the boundary of both of the Poincare models. Any two points 
on the surface are connected by a unique geodesic. The unique geodesic that connects these 
fixed points is called the boost's invariant geodesic. 

2. Rotations : these isometries are the analogues of Euclidean rotations. They have one fixed 
point which does not lie at infinity. 

3. Inversions : reflections across the invariant geodesies of boosts. 

4. Odd boosts : an even boost followed by an inversion through it's invariant geodesic. 



One of the important facts about the pseudo-sphere is that the isometries listed above, are given by 
Mobius transformations in the Poincare disk and half-plane: 

(' = g{() = ad — be = 1, a,b,c,d gH , (3-3) 

in the case of even boosts and rotations (in the half-plane), and 

C' = 5(C) = ^7^ ad-bc = -l, a,b,c,dGli, (3.4) 
cC -I- d 

for inversions and odd boosts (in the half-plane). This means that an isometry that is obtained by 
a sequence of isometries, is the product of the appropriate matrices. 

Many statements about the invariant geodesies of even and odd boosts are going to be made in the 
following discussion. They are usually easier to consider in the half-plane model, because in it the 
positive real axis is a geodesic, therefore it is always possible to chose coordinates such that it is 
the invariant geodesic of the boost under consideration. In this frame of reference the matrix that 
corresponds to a boost is diagonal: 

-|- even boost , , 

- odd boost ^ ■ ^ 




As the curvature is negative, the flow upon the pseudo-sphere is hyperbolic. This means that adjacent 
trajectories diverge exponentially fast from each other. In spite of this, motion on the whole surface 
is integrable |20|. According to ||2^, hyperbolic flow on a compact surface is chaotic. Therefore, 
motion in any compact surface on the pseudo-sphere is assured to be chaotic. One such choice are 
tiling billiards. Such billiards are the primitive cells of lattices of geodesies on the pseudo-sphere, 
as can be seen in Fig. |l[ The word "billiard" refers to an enclosure on the pseudo-sphere with any 
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boundary conditions. In this paper two cases are going to be concentrate upon, namely periodic and 
Dirichlet boundary conditions. 

The situation presented above is very different from the one arising in the Euclidean plane. In that 
case trajectories diverge only linearly on the whole surface, so that in order to ensure the chaotic 
nature of the motion complicated boundaries have to be chosen to introduce the necessary exponential 
divergence. 

Having a billiard that is a primitive cell of a lattice of geodesies, means that there is an underlaying 
discrete group of isometries, whose fundamental domain is the primitive cell. A fundamental domain 
of such a discrete group is a domain, that under the action of all of the different members of the 
group tiles the whole surface perfectly. The simplest group consists of only even boosts (in addition 
to the identity), the fundamental domain of which is a polygon with periodic boundary conditions. 
If the fundamental domain is to be a billiard with hard walls, than the group consists of even boosts, 
odd boosts, reflections and rotations. 

An interesting feature of tiling billiards on the pseudo-sphere is that their periodic orbits are closely 
related to the boosts of the underlying group. This situation also exists in tiling billiards on the Eu- 
clidean plane (with translations instead of boosts). The important difference between these systems 
is that in the latter case the underlying group is Abelian, whereas in the former it isn't. In addition, 
on the Euclidean plane there are continuous families of periodic orbits, whereas on the pseudo-sphere 
all periodic orbits are isolated. 

Consider free motion inside the fundamental domain, the billiard. Since the motion is free until the 
boundary is hit, it is on geodesies. Once one of the edges is hit the trajectory proceeds according 
to the boundary conditions. Again, the simplest case is that of periodic boundary conditions, in 
which the trajectory re-enters from the edge identified with the one it hit, on a new geodesic. A 
more mathematically complicated case is the case of reflections, in which the trajectory is reflected 
specularly and continues on the appropriate geodesic. 

There is another way to describe the motion, in which the fact that all the different copies of the 
billiard, under the discrete group of isometries, tile the whole surface. Beginning on a geodesic 
segment in the billiard and travelling on it with constant velocity, eventually one of the edges is hit. 
Then, instead of staying in the billiard on a different geodesic, the original geodesic is continued into 
the billiard copy across the edge. If the boundary conditions were periodic this copy is identical to 
the original billiard, and if the edge is a hard wall, the copy's orientation is reversed with respect to 
the original billiard, since it is obtained by reflection. 

Now, the next segment of the trajectory is on the original geodesic, but in a different copy of 
the billiard. All this process will have to be repeated once another edge is hit. Continuing the 
construction ad infinitum, the trajectory is described by the collection of the different copies of the 
billiard through which the original geodesic passes. 
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Recapitulating, instead of describing the classical trajectory in the usual way, as a collection of 
geodesic segments, all of which lay inside a limited area of the surface, it can be described with one 
geodesic that enters and exits copies of the original billiard. 



The above procedure is standard practice in polygonal billiards (see for example ||3^). The next 
step's presentation is due to Bogomolny et al. [31|: 



1. Consider one of the geodesies described above, that enters and exits many tiles of the tessel- 
lation; 

2. Each tile along it is labeled by the group isometry with which it was obtained from the funda- 
mental domain: {I, gi,g2,P, • • •} {I is the first segment, the one in the original copy); 

3. To reconstruct the "physical" trajectory, the geodesic must be folded back into the original 
fundamental domain. This is done by copying each geodesic segment by the isometry inverse 
to the one that created the copy it is in. These isometries are: |/, Oi^ , :P^^ j • • - j; 

4. After this operation most of the geodesic segments will be dispersed around in a complicated 
manner. But, for special trajectories, one of the segments (the one that was in the copy that 
was suggestively labeled p) will match the first segment exactly in location and in orientation. 
In this way: 

an isometry that copies a geodesic segment onto a segment with identical orientation 
of the same geodesic has been found. Therefore: 

(a) The geodesic segments constitute a periodic orbit; 

(b) This geodesic is the invariant geodesic of p; 

(c) On the pseudo-sphere p has to be a boost, for it has an invariant geodesic. (It 
can't be a pure reflection because the orientation is the same.) 



The arguments given here are applicable to tiling billiards on the plane too, where the role of the 
boosts is taken by simple translations. The qualitative difference is that on the plane, periodic orbits 
are not isolated, whereas on the pseudo-sphere they are, due to the exponential divergence of close 
trajectories. 

The distance a point is transfered by a boost is closely related to the it's distance from the boost's 



invariant geodesic. If is an even boost p4|| 



1 1 

sinh-dg(C) = cosh(i^(C) smh-L„, (3.6a) 

2 2 

where dg{C) = df-^g(^Q is the distance of a point from it's image under g, and (C) is the distance of 
the point C from the invariant geodesic of g measured along a perpendicular geodesic (see Fig. 0). 
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A similar relation can easily be obtained for odd boosts: 

cosh -(ig(C) = cosh cosh -Lg. (3.6b) 

3.1 Green's function for fundamental domain billiards 

The Green function is calculated exactly by the method of images, from the free Green function, 
by virtue of the underlying discrete group of isometries. In half-plane coordinates, Schrodinger's 



equation reads [gO| 



_^/f^ + ^U = E*. (3.7) 

Distance is measured in units of R, and from this point on, energy will be measured in units of 2mK^ 
so that = if m = ^ then Ti = \ . In these units, it turns out that the free retarded Green function is 



given by [g^ 



G^(C,C') = ^Q-i-.fc[coshdc,d , (3.8) 
where E = k"^ + \ and Qi[x\ is Legendre's function of the second kind. 

The billiards under consideration are the fundamental domains of groups of isometries. Therefore, 
all the copies of the billiard tile the surface perfectly. In addition, as all the copies are identical, 
a solution of Schrodinger's equation in one, solves it in all the other tiles, but for a coordinate 
transformation. This means that Green's function has to be invariant under the action of the group 
isometries {gi € F). A simple way to obtain such a function is to use the method of images and free 
Green's function: 

GE{Q.C') = Y.XgG%{C.9{C')) ■ (3.9) 

where Xg = (±1)^""*^ °f ^ (+1 is the choice for Neumann and periodic boundary conditions and 
— 1 for the Dirichlet case). 

Now the exact retarded Green function for fundamental domain billiards can be written using (| 

G'i?(C,0 = ^Ex9^?-i-.fch^^^C,^;(C')] ■ (3-10) 

Eq. 3.10 has an intuitive explanation. It gives the Green function between Q and Q' as a sum over all 
the classical orbits between them. To see this one should describe a complicated orbit between these 
two points using the scheme described already in the previous subsection. The distance between C, 
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and one of the copies ^(C') is nneasured along the unique geodesic that connects them. As ^(C') 
is completely equivalent to this geodesic is an (unfolded) classical orbit between C, and This 
is fundamentally different from other quantum mechanical problems, where all Feynman paths are 
required for the calculation of Green's function. This property, that holds when tiling is possible, 
makes the semiclassical approximation exact for the problems studied here. 

The sum over the group members of (3.10) can be arranged in a suggestive manner. Any group 
can be divided into conjugacy classes |^. A conjugacy class is defined as the collection of all 
different group members that are similar to one another. All the members of a conjugacy class can 
be generated from one representative via: 

g' = h-^gh h G T/Tg , (3.11) 

where T/Tg is the sub-group of T, whose members are unity and the members that do not commute 
with g. The reason for excluding Tg (the sub-group of the group members that commute the g) is 
that each member of a conjugacy class should be accounted for only once. 

It may be shown that there is a one-to-one correspondence between the conjugacy classes of boosts 



and the periodic orbits [g5|, This means that the list of boost conjugacy classes is the list of 



periodic orbits. Therefore, it has the form: 

{...,pI,...} , 

where r = 1,2,3, ... are repetitions and {pi} is the list of primitive periodic orbits of the billiard. 
Note that p and p^^ can only be conjugate if there are reflections in the group, i.e if the periodic 
orbit is self retracing. 

The exact retarded Green function is given by (3.10). Rearranging the sum over boosts in complete 
boost conjugacy classes, one obtains: 



Ge{C,C') = i^T.X9Q 



2tt ^ --^y-'-a" 

a 

not a boost 



ik 



cosh d(^^g((;') 



+ 



^ CO 

{p} r=iher/r^ 

where Z^^^i is the sum over all the different boost conjugacy classes, and J^her/Fg, i^ the sum 
over the different members of each class. 

For a billiard with periodic boundary conditions, the group of isometries, whose fundamental domain 
is the billiard, contains only the identity and even boosts. Retarded Green's function is given, in this 
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case, by: 




cosh d, 



'C,ft-V/i(C') 



(3.13) 



where J2{p} E^i 



is the sum over all the different even boost conjugacy classes. 



For a billiard with hard walls, the group whose fundamental domain it is, is generated by reflections 
at the edges. This means that it will include the following types of isometries: 

1. The identity /. 

2. Inversions {1} - reflections through the edges of all the tiles of the tessellation. 

3. Rotations {TZ} - rotations around all of the vertices of all of the tiles in the tessellation. 

4. Singular boosts - hyperbolic isometries that commute with some reflection. They correspond 
to periodic orbits that run along edges of the billiard. These periodic orbits are invariant to 
reflections across those edges, so singular boosts always come in pairs of odd and even. 

5. Odd boosts - correspond to periodic orbits that hit an odd number of edges of the billiard. 

6. Even boosts - correspond to periodic orbits that hit an even number of edges of the billiard. 

As before, the boosts can be divided into conjugacy classes, whose list is the list of periodic orbits, 
which are repetitions of primitive periodic orbits. Since there are inversions in the group, there are 
two types of primitive periodic orbits: those that bounce an even number of times and those that 
bounce an odd number of times, during one period. The former will be denoted by {pi}, and the 
latter by {qi}. 

It is convenient to separate the singular boosts from the rest. The previous argumentation about 
hyperbolic isometries applies to them, so they too can be ordered according to conjugacy classes. 

One of the peculiarities about the orbits that correspond to their conjugacy classes, is that for every 
even primitive singular orbit, there is an odd one. This is why these conjugacy classes can be ordered 
in pairs: {s} and {X^s}, where Is is an inversion that commutes with s. 
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Now (3.12) can used again to write the exact retarded Green function of the problem: 



2tt 



cosh d, 



CMC) 



+ ^EQ- 



-ik 



cosh d, 



+ 



EE E Q 

{p} r=i her/Fp 



5 *^ 



cosh d, 



'(:,h--^p-hi(') 



+ EE^^E Q-^- 



ik 



{q}r=i h&r/r. 



cosh d, 



+ 



OO 

+EE E [Q-^,-.k 

{s} r=i her/r. 



cosh d, 



'C,h-is'-/i.(C') 



+ e<3_i_jfe coshd^^;,-ij^^r;,(^,) 



(3.14) 



where e = — 1 for Dirichlet boundary conditions, and e = +1 for Neumann boundary conditions. In 
this work only the former will be studied, because the generalization for the latter is straightforward. 
The term "hard wall billiard" will be used from this point forth to describe a "Dirichlet boundary 
conditions' billiard" . 



3.2 The density of states and Selberg's trace formula 



The density of states can be calculated exactly by tracing Green's function. For simplicity, a billiard 
with periodic boundary conditions will be considered initially. Tracing (3.13) yields the density of 
states that takes the form [pO| : 

A i L °° pikrLp 

+ = ta„h,. + « 1:^1: } . (3.15) 

Up} '■=1 2 



where A is the area of the billiard. Eq. 3.15 is the Selberg trace formula [g5|. 

Using Eq. 3.15, one can derive Weyl's asymptotic expansion for the counting function of eigenvalues 
(BV VII.36): 

N{E)^^(e-1) as ^^00 (3.16) 



4vr 

It turns out that all other terms in the series are exponentially small as ^ 00. 

The computation of the density of states in hard wall billiards is very similar in principle, so it will 
not be presented. A detailed computation appears in [|^. The result of tracing (3.14) is: 

p{E)=p{E) + ppo{E) , (3.17) 



16 



where (BV M.U): 

{A C 
— tanh nk — — - + 
47r 8-Kk 

1 v-^ 1 dxcoskx ( , m,rX , 2;\ , , 

,5, ^hS|- — - 2) • 

in which {;^} denotes the list of the billiard's corner angles, £ is the perimeter of the billiard, and 

where {p} denotes the list of all primitive periodic orbits, even, odd or singular. In the last equation, 
the definitions: 

Lp even p.p.o [ even p.p.o 

Up = Lq + in odd p.p.o 7p = S odd p.p.o (3.20) 

2Ls singular p.p.o , I Ls/2i singular p.p.o . 

have been used. As for billiards with periodic boundary conditions, one can derive Weyl's asymptotic 
expansion for the counting function for eigenvalues for hard wall billiards using Selberg's trace formula 
(|3 17| ). For Dirichlet boundary conditions the result is (BV VII. 50): 

N(E) — E - - — + y —(mr -—]+■■■ asE^oo. (3.21) 

^ ^ 47r 47r 127r r4^, 24 V rrir J 



3.3 Spectral determinant 

Because the Selberg trace formula, that holds for billiards on surfaces of constant negative curvature, 
is exact, while Gutzwiller's trace formula is valid only in the semiclassical approximation for hard chaos 
systems, many of the semiclassical expressions cited in Section |2| are exact for tiling billiards on the 



pseudo-sphere. For example, Eq. 2.8 is an exact equality for the spectral determinant if one uses the 



definitions in (3.20). Therefore, ( |2.9| ) is also an exact equality. On the other hand, resumming the 



problematic sum over primitive periodic orbits in ( P^ ) yields only a semiclassical approximation to 



the exact result. Eq. |2.11| must be adapted to the unit convention and billiards used here. In the 
units being used h = 1 so S = kL, where L is the length of a trajectory, and the wavenumber k 
plays the role of ^, becoming large in the semiclassical approximation. 
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Replacements for Eq. 2.12a and Eq. 2.12b are also required: 



Ak 

i^{n,E)^i^{k)^L^- — , (3.22a) 

Q\k\)^^^K^^.^. (3.22b) 

where = Jp^p- The result of the resummation for a particular hard wall billiard (defined in 

Fig. |l|) is shown in Fig. |3|, together the exact zeros of the spectral determinant, calculated numerically. 
One can see that the agreement is excellent even at the lower end of the spectrum. 



4 Green's function and eigenfunctions on the pseudo-sphere 



In this section the diagonal part of Green's function for billiards on a surface of constant negative 
curvature is calculated. More precisely, it is shown that for C = C' the general results (3.9, 3.10, 3.12 
and 3.14) can be written in the form (1.1), namely, as a numerator that depends on the coordinate ( 
over a denominator that is the spectral determinant. Both numerator and denominator are expressed 
in terms of contributions of various classical orbits. In this way, a comparison to standard semiclassical 
results is made possible. The eigenfunctions can be calculated directly from the expression of Green's 
function that is obtained in this section, with the help of (|2.2|). 



' m 



At this point a clarifying remark is due, because Green's functions are ill-defined for r 
Green's functions satisfy the equation: 

(A + ^) Gs(r,r') = <5(r-r') , 

and in two dimensions this necessarily means that the Green function has to have a logarithmic 
singularity as r ^ r'. This singularity can be removed by subtracting the divergent part: 



lim 
laHO 



1 



Gsir, r + a) - — In dr,r+a 
zvr 



GE{r,r) 

This subtraction does not affect the spectrum, since it does not depend on k. 



As in Section |3.2[ it is instructive to start by considering a billiard with periodic boundary conditions 
on the pseudo-sphere. A generalization for the case of a billiard with Dirichlet boundary conditions 
is presented in Section Beginning with the exact (retarded) Green function for ( = given by 
(3.13): 



GEiCC) 



-1 
2^ 



-1 

2^ 



EE E Q- 

p r=iher/rp 



-ik 



cosh 



(C) 



(4.1) 
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Notice that had the above mentioned definition of G£;(r, r') not been used, the right hand side of 
the last equation would include a logarithmic singularity, because (BV G.20): 

i„j^ [coshd] ~ — Ind as d ^ . 

In Section ^ it was noted that the sum over hyperbolic isometries in (4.1) contains information about 
the lengths of the periodic orbits of the billiard. In this chapter more information is going to be pried 
out of it, in order to sum the repetitions. This will be achieved with the aid of ( 3.6a| ). This simple 



relation from hyperbolic geometry will facilitate the separation of the two different contributions to 
the distance between a point, and it's image under a hyperbolic isometry, g{C)- One contribution 
to this distance is the distance of C from the invariant geodesic of g, (C)- The second contribution 
is the length of the periodic orbit that corresponds to the isometry, Lg. The treatment of this latter 
contribution is going to be of vital importance in the ensuing discussion. 



4.1 An exact formula for eigenfunctions 



In order to obtain a transparent formula for Green's function, Legendre functions, which constitute 
it, are expanded in powers of e^^f (see Appendix"^ ): 



where 



Q_i_,fc[cOshd,(C)] =Y.^rn [k,dj{C) 
m=0 



g-(i+m-ifc)L9 



(4.2) 



X dm— 2n 



n=0 ^ ' 
1 

cosh^ d^iC) 



drj 



1 ., 1 

n H ik, — 

2 '2 



^cosh^4(C), 



in which B{x,y) = ^^^(J^+y) ^ets^ function, and a„i [x,y\ is the product two real polynomials 

of order m, one of x and one of y. 



(4.3) 



Returning to (4.1) and using (4.2), one obtains: 



GE{CX) = —Q-l-rk[^\ + 



27r 2- 



{p} r=l heV/Tp m=0 



(4.4) 



19 



where was defined in (4.3). 



It was already stated (in the paragraphs before Eq. 3.12) that all the members of a conjugacy class are 
related to one periodic orbit. That fact, together with the explanation in the paragraphs preceding 



(|3.6a|) , leads to: 



L 



rLr- 



(4.5) 



Consider the set of distances < ci^_i r/j(C) f of the point C from various invariant geodesies. 

The most obvious point to notice is that. 



because p"^ and p have the same invariant geodesic. 



(4.6) 



Another point worth mentioning is that the set of distances < dj^_i ^(C) \ is, actually, the set 

of distances of ( from all different images of the invariant geodesic of p§^. This is true because if 
two boosts {gi Sz §2) have the same invariant geodesic they must be repetitions of one another. The 
reason is the one-to-one correspondence of boosts and periodic orbits: if two periodic orbits overlap, 
they must be repetitions of one another, because different classical trajectories never cross in phase 
space. In a billiard with periodic boundary conditions this means that gi and §2 are identical§^. 

Now the summation over the distances from the invariant geodesies of members of a hyperbolic 
conjugacy class may be given a geometrical interpretation. What one actually does in this summation, 
is to take into account the distance of a point from all of the images of the invariant geodesies of 
the class representative. 



The next aim is to derive Eq. 4.1C, 4.11 and 4.12, which express the exact Green function in terms 



of sums over periodic orbit lengths. This will be achieved by using Eq. 4.5 and 4.6 in Eq. 4.4 



GEiC,C) = ^Q^l-^k [1] + 



+ 



2lT 



{p} r=l h&r/Fp m=0 



(C) 



3— (^+m— ifc)rLp 



(4.7) 



Note that the contribution of a single periodic orbit (p) has been decoupled into two separate parts - 
the perpendicular distance, dj^-iphiO' the length of the periodic orbit, Lp. To achieve this in a 



^^The only exception to the rule occurs in hard wall billiards because they have self retracing periodic orbits. 
Then the list of distances is doubled, because images of are in the conjugacy class. 

§^In the Dirichlet case gi = if g corresponds to a self retracing periodic orbit 
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typical chaotic system one usually has to resort to the semiclassical approximation, while here it is 
exact. 



Exchanging the order of summation, a geometric series is encountered: 

^ ^-{\+m~ik)Lp 



+m—ik)Lp 



r=l 



(4.8) 



This enables one to write the retarded Green function (4.7) in the form: 



— 1 —1 °° r 

Ge{Q.Q) = ^Q-i-.. W + E E •^'n \k,dUr>h{0\ S^^\L,,k) , (4.9) 



27r 



{p} her/Fp m=0 



that is given in terms of the lengths of the primitive periodic orbits of the system. Taking a common 
denominator and multiplying numerator and denominator by e~*'^^(^\ one obtains a convenient 
exact expression for the retarded Green function: 



GEiCO 



AiE) 



where 



-1 



^iE)^—Q_r_^,il)A{E) 



(4.10) 



(4.11) 



in which A{E) is the spectral determinant, as it is given exactly by ( ^ ) using the definitions of 
(3.20) and the replacements defined in Section 3J in order to utilize Berry and Keating's results. In 
addition: 

M,oic,E)^—Y: E E-^-[^'4-WC)] A(p'-)(ii;) . (4.12) 

{p} her/Fp m.=o 

A(P>'")(£;) is defined in (2.17). According to Eq. ^ and |2l8| , as explained in Section U, A{E) 
and A(^''"^)(£') are given exactly by sums over composite orbits: 



Cfj, and cj^'™'^ are given by their definitions Eq. 2.10b and 2.19b , with 3.20. It turns out that all 
of these coefficients are real numbers (see [Appendix CI ). The L^'s were defined after (3.22a), and 

^ti,P — + Lp. 

The criterion for the separation of the numerator in Eq. |4.10| into two parts, is that tracing the first 
term gives the smooth term in the density of states to (3.15), and tracing the second one gives the 
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periodic orbit contribution to (3.15). 



There are three sorts of contributions to Eq. 4.10. The first is from short non-periodic orbits, and it 



is given by J\f{E), as it results from the identity. The second and third contributions are related to 
periodic orbits and are given by J\fpo{C,E). The first contribution to this function is related to the 
actual periodic orbits in the billiard, while the second is related to their images. The reason is that 
a "real" periodic orbit is made up from only a few segments of invariant geodesies of a conjugacy 
class's elements. The invariant geodesies of the rest of the conjugacy class make up "images" of 
the periodic orbit, as they never intersect with the original billiard. "Real" periodic orbit segments 
generally have relatively small dj^^-ip^iO' while the images of them tend to have larger d-^^i^^^iC,). 



Eq. is a manifestation of the results of Fredholm's theory, as presented in flj]. Unlike in 



the cases considered there, the exact Green function is written as the ratio between two sums over 
periodic orbit lengths without the need to resort to a semiclassical approximation. The precise nature 
of the connection with Fredholm's theory has not been worked out yet for this case. Once this issue 
is resolved, Fredholm's theory will supply a prescription as to the correct reordering of terms in the 
series, so that they converge absolutely to the exact result. 



The eigenfunctions densities can be found from (4.10) with the help of (p^). The poles of Green's 



function are the eigenvalues of the Hamiltonian, and the exact eigenfunction densities are the residues: 

\M0\^ = Res \ / = \ / , (4.13) 

where A'(^Ea^ denotes a derivative with respect to energy of the spectral determinant, evaluated at 
-Sa- 
in order to make sense of the infinite sums over the periodic orbit lengths resummation is required, 
as they do not converge absolutely. This is performed in the Section |^. 



4.2 Hard wall billiards 



A formula for Green's function for tiling billiards with Dirichlet boundary conditions on the pseudo- 
sphere that is analogous to Eq. 4.10| - 4.12 will be developed in this subsection. The differences 
from the periodic boundary conditions' case arise from the fact that now there are several new kinds 
of transformations: inversions (X), rotations [TZ) and odd boosts. In principle, the derivation of a 
formula for Green's function is exactly the same as the previous derivation. As before, the exact 
Green function can be written as the ratio between two functions that involve summation over all 
the periodic orbits of the system. The differences are of technical nature only. 



As before, one may separate Green's function into part that is related to the periodic orbits and to 
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a part that is not (see Eq. 3.12). The criterion for the separation is the same as mentioned after 
(|4.12| ), namely tracing the first term gives the smooth term in the density of states (3.18), and 
tracing the second one gives the periodic orbit contribution (3.19). According to (3.14), for the 
Dirichlet boundary conditions case, the latter is given by: 



-1 f °° r 



\, {p} r=i her/Vp 



+ 



{q} r=i her/r, 

oo 

+EE E 

{s} r=l her/Fs 



+ 



+EE E (-ir^?-H^ 

The three types of boosts that exist have been separated: 



(4.14) 



1- J2{p} Z^J^i is a sum over the classes conjugate to powers of primitive even boosts. 

2- J2{q} S^^i is a sum over the classes conjugate to powers of primitive odd boosts. The factor 
(—1)^ is necessary, because ifr is odd/even hq^h~^ involves an odd /even number of inversions. 

3- Z]{s} is a sum over the classes conjugate to powers of primitive singular boosts. These 
boosts are the boosts whose invariant geodesic lies along an inversion line of the group, which 
are images of the boundary of the billiard (see BV Appendix M). This means that [s,Ts] = 
[TsS,Ts] = 0. This is why both s and T^s must be included, and why the contributions have 
opposite signs. 



For the primitive even boosts a contribution identical to (4.2) is obtained, while for the primitive odd 
boosts the result is (see |Appendix~B| ): 



(-l)''Q_i_ifc [cosh4gr;,-i(C) = 

oo 

= (-ir E ((-i)T-^" 



- (l+m— jfc)rL5 



(4.15) 



m=0 



where J^„ 



A;,d^_i^^(C) was defined in Eq. 4.2. For the primitive singular boosts a slightly different 



expression is obtained (see Appendix B): 

Q-l-ik COsh4sr/i-l(C) COsh4jsr;j-l(C) 



m=0 



3-(|+m-jfc)rLs 



(4.16) 
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The important fact is that, as in (4.2), the contribution of the periodic orbits to the exact Green 
function is now expressed as a sum over periodic orbits in a form that admits summation over 
repetitions. Since Eq. 4.2, 4.15 and 4.16 have a similar form, the same steps that were undertaken in 
the previous case apply here. After summing over the repetitions and taking a common denominator, 
the expression ( |4.10D is found, but J\f{E) and Afpo{C, E) are of a different form. The function Af{E) 
will be treated later, while: 



-1 



Ko{c,e)^—Y: E E^^h4-v/.(0jA(^'™)(i?), 

{p} her/r^ m=o 

where J2{p} denotes summation over all primitive periodic orbits, even, odd or singular, and: 



(4.17) 




P / s 
p = s 



Using the definitions of (3.20), A{E) is given by (U) and A(P''")(S) by (|]l|) 



One can see that the additions to ( |4.12| ) required for the calculation of the contribution of the 
primitive periodic orbits to the eigenfunction density in billiards with hard walls, do not change it's 
nature, adding only contributions of other types of periodic orbits, namely odd and singular ones, 
and their corresponding amplitudes and phases. 

The contribution of non-periodic orbits is given by the contribution of the non-boost isometries 
of the group, whose fundamental domain is the billiard under study. In the case of billiards with 
periodic boundary conditions there was only one such isometry, namely the identity, that corresponds 
to zero length orbits. In the case of billiards with Dirichlet boundary conditions there are many such 
isometries besides the identity, all of them rotations and inversions, and they correspond to orbits 
closed in real space, not in phase space. According to Eq. 3.14, the contribution of these isometries 
to Green's function is: 



G£;(C,C) = ^ 



E Q-i-.fc[coshdx(C)] + 

ions 

+ E Q-l-^k[^oshdniO]\ 

Rotations J 



Inversion, 

I 



(4.18) 



■R 
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This means that M (E) of ( |4.11| ) should be replaced by: 



2tt 



Inversions 

I 



Q_i_.Jcoshdx(C)] + 



+ '3-i-ifc[coshd7e(C)] 

Rotations 

n 



(4.19) 



The main difference between ( |4.11|) and this result, is that the former is a constant function of 
position, whereas the latter is not. 



5 Resummation 



As it stands, Eq. 4.1 and all equations derived from it are not very useful because of the same 
convergence difficulties (|2.6| ) presented. The reason is the exponential proliferation of the number 
of images of (" with the distance dg{(). When large enough, the distance d = dg{() is approximately 
the distance of the the g'th tile from the billiard. The number of tiles at such a distance is roughly 
proportional to the circumference of a circle of radius d around the original billiard. On the pseudo- 
sphere this circumference is proportional to e'^ as d ^ oo. On the other hand, the behavior of 
Q_i_,;fc[cosh(i] for large d is: 



Q_i_jfc[cosh(i] ~ e(*'' z)"^ as ci ^ oo , 
so the net contribution of images at a distance d to (4.1) is proportional to 

e('''+^2)d as d ^ oo . 



(5.20) 



(5.21) 



This means that the sum in (4.1) is at best only conditionally convergent, which is why resummation 
is needed to make use of it. 

Notice that most of the exponential proliferation that causes the convergence problems comes from 
the exponential proliferation of classes conjugate to primitive isometries, which is equivalent to the 
exponential proliferation of long primitive periodic orbits. This is because the number of periodic 
orbits of length d is [|4[^: 



N{d) 



d 



as d 



oo 



(5.22) 



which means that the number of repetitions of short primitive periodic orbits is sub-dominant and 
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that the number of terms in the sum over a conjugacy class's elements grows sub-exponentially. 
Resummation following Berry and Keating |^ will be applied to get a meaningful result from this 
sum. It will result in an approximation on the level of the semiclassical approximation 



Due to Eq. |5.21| , the sum over {p} in Eq. |4.1C1| does not converge absolutely for real k's. On the 
other hand, if: 

9=A; > ^ (5.23) 

it does. It will be assumed that J\f(^Ci E{k)j is analytic in a sufficiently wide strip around the 
axis to resum the problematic sum by analytic continuation from the region described by (5.23) to 
the region of physical interest Qk = 0. 

In the previous section, a formula in terms of periodic orbits was obtained for the exact expression; 



Gi^(c,c) = E 



E{k) - 



A(^E{k) 



(5.24) 



Resummation of the denominator was discussed in Section [3.3[ so only the numerator will be treated 
explicitly in this section. 



Since E{k) = /c^ + 4 = E{—k), and since 



A(E{k)) = i^(E{-k) 



then: 



M((:,E{k)) =M(Q,E{-k) 
This exact symmetry is crucial for the analytic continuation method of resummation to work 

Using Cauchy's theorem one may write: 



(5.25) 
(5.26) 



M(c.E{k) 



2-111 Jc z 



''-^{z,k)N(c,E{k + z)) . 
Leading to: 

-cx)+i(l/2+e) 

AA(c, E{k)) = ^ j ^7(z, k) [aA(c, E{z + k)) + AA(c, E{z - k 

oo+j(l/2+e) 

where the ad hoc choice of Berry and Keating |^ will be used: 

7(z, A;) = e ■^^'^ . 



(5.27) 



(5.28) 



(5.29) 
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As shown in |/\ppendix~D , Eq. 5.28's leading order approximation is given by 



where 



M[E{k)) ~ {A{E) + iAi{E)) | as A; ^ oo . 



A{E), Ai{E) G M are given by: 

A{E) + iAi{E) = ^ C'^e-*'^^(-^)+^^^'' Erfc 

where ^^(fc) and Q'^{K,k) are given by (3.22a) and (3.22b). 
In addition: 



k 



2Q^iK,k) 



Ko{CE{k))^^^Y. E A(^''')(C,^)1 as A: 



CXD 



where: 



in 



A(f''^)(C, A:) = ^ Ajf''^)(C, fc)e-*"^(^)+*'=^''-'' Erfc 
which i^^p{k,C) = ^M,p + l^cosh^(i^p^_i(C) - ^A;^^, and: 



2Q^iK,k) 




E 



sinhd 



hp3+^h- 



(C) 



as A; ^ cx) 



where C^^'^'^ was defined in (|2.23|), and the definitions of (3.20) for even boosts were used. 



(5.30a) 
(5.30b) 



(5.30c) 



(5.31a) 



(5.31b) 



(5.31c) 



Since both sides of ( |5.24j) have only simple poles at the eigenvalues of the Hamiltonian, the residue 
at E = Ea is |V'a(C)| ■ Hence one obtains: 



iV'a(c)r ~ 1^^.(01 +\M0\;o as^^oo 



where 



\MC)\ 



■2_ M{C,Ea 

A'{E^) 



(5.32) 



(5.33) 
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(5.34) 



Notice that on the spectrum; 

M{C,Ea) = ^Ai{Ea)tanhTrka ■ (5.35) 

In (5.31b) an effective truncation of the problematic sums over periodic orbits (=conjugacy classes) 
in ( |4.12| ) was obtained. The cut-off is centered around 

Cip{kaX) = 0, (5.36) 

that corresponds to 

L^,, + lncosh2 4^,_.(C) = ^ , (5.37) 

which is equal to the distance traversed during half the Heisenberg time. There are two sorts of 
contributions to the sum over periodic orbits. The first is from geodesies that pass through the 
billiard. Segments of these geodesies constitute the "real" periodic orbits of the system, so they 
are characterized by relatively small c?^p^-i's. The second contribution is from geodesies that never 
intersect the billiard. They only pass through copies of it, so they are characterized by relatively large 
values of c?^p/j-i's. This is why the first contribution can be associated with periodic orbits, and the 
second with their images. These two sorts of contributions are translated to two sorts of cut-offs. 
The first is the cut-off of long periodic orbits, which is encountered in typical chaotic systems too, 
while the second is the cut-off of distant images of short periodic orbits. 



5.1 Resummation for hard wall billiards 



As for billiards with periodic boundary conditions, in order for ( 4.10 ) to be meaningful for hard wa 



billiards, resummation is required, because Eq. 4.17 and 4.19 involve sums that do not converge 
absolutely. The same scheme Berry and Keating Q introduced may be used again. The only quali- 
tative difference results from Weyl's asymptotic expansion for the counting function of eigenenergies 
(|3.21j ) because it's high order coefficients are not identically zero, as they were in ( [3.16|) . Despite of 
this, in the leading order approximation, there is no difference. The result of the resummation is, to 
leading order: 



I {p} /ler/rp J 



as k 



oo 



(5.38) 
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where: 



Eq. 5.39 includes the contribution of Eq. 5.31c (even boosts), odd boosts: 



3q 



j=0 



sinhdhqj+ih-^ (C) 



as well as singular boosts: 



j=0 



sinhd 



(C)J 



as A; — > cx). 



where the coefficients C^'''"''' and C;^*'"''' are those of (p.23|) with the definitions of (3.20). 



(5.39) 



(C)-{i+i)i, 

AL'''''^(C,fc) ~ J-Vy(-l)^'+^Ci'''^')e*(^'+i)^'^ =^ ask^ cx),(5.40a) 



(5.40b) 



As to the contribution of the non-periodic orbits, (4.19) leads, after the analytic continuation, to the 
expression that should replace (5.35) : 



Ai{Ea)tanhTTka , 

A/ (.-c/Q, <; j = h 



(5.41) 



E Q- 



[cosh dj(C)] Erfc 



Inversions 

I 



+ 



+ E Q i_,fcJcoshd7e(C)]Erfc 



Rotations 

n 



2Q\K, k. 



6 Numerical calculation of eigenfunctions 



Eq. 4.1C for billiards with periodic boundary conditions (with 4.11 and 4.12| ), or for hard wall billiards 
(with 4.17 and 4.19), is exact. It cannot be used, however, for the calculation of eigenfunctions, since 
both the numerator and the denominator are expressed as series that are not absolutely convergent. 
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In what follows, resummed expressions obtained from the above equations, namely Eq. 15.321 with 5.38 



and 5.41 will be calculated and compared with exact eigenfunctions found numerically. In particular, 
an arbitrary slice through the two lowest eigenfunctions of the hard wall billiard defined in Fig. |l| is 
calculated. This billiard tiles the pseudo-sphere perfectly by reflections, because all of it's angles are 
of integer fractions of vr At such low energies, the semiclassical approximations used in the 

resummation are not expected to be of very high accuracy. The reason higher eigenfunctions were not 
investigated is that a large number of isometries has to be used in the calculation. There are billiards 
for which it is very easy to obtain large numbers of isometries. These billiards are called "Arithmetic 
Billiards", and the reason for not choosing one of them is that some of their physical properties, such 
as the energy spectrum and the spectrum of classical periodic orbit lengths, is non-generic pi 



The results are shown in Fig. |5|and ^. Presented are the exact cut through the eigenfunction density, 
the contribution of the the non-boost isometries of (5.41), and the corrections to this due to the 
periodic orbits of (5.38). The contribution of the non-boost isometries accounts for most of the 
features of the eigenfunction density, while the addition of the periodic orbits provides only a small 
correction. 

One can see that even after the inclusion of all of the isometries the sums truncated by resumma- 
tion require, the eigenfunctions are still not reproduced perfectly. A measure of the quality of an 
approximation of the function f{x) by fapprox{x) is given by the average squares deviation: 

/ DX — I '^•^ ~ fapprox{x)\ , . s 

J dx\f{x)\ 

where in this case f{x) = |^'(C)|, and the integral is taken along the cut (see Fig. It is inter- 
esting to see the difference in the quality of the approximation without including the periodic orbit 
contribution, (Di), and including it, (-D2)- For the ground-state (see Fig. |5|): 

(Di) « 0.023 , 
{D2) « 0.015 . 

For the first excited state the result is: 



{Di) ^ 0.043 
{D2) « 0.030. 

These numbers show that the quality is fairly good even without including the periodic orbits' 
contribution. Including this last contribution improves the approximation appreciably. The reason 
the fit is not perfect is that at such low energies the number of periodic orbits that contribute to 
the truncated sums is very small. The cut-off region extends all the way from Lp = to Lp = 2.00 
for the ground-state and to Lp = 2.46 for the first excited state. This means that almost all 
the periodic orbits that contribute, in both cases, to the eigenfunction density are from the cut-off 
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region. Therefore, one cannot expect that the truncated sums in (5.39) and (5.41) should do too 
well, because in deriving them, terms of the same order were neglected, namely terms supported by 
the cut-off region exclusively. For such low levels it is essential to include high order corrections to 
the truncated sums, and although it is possible to do so in an obvious manner, it was not done to 
date. 



7 A comparison with previous results 



One observes that the structure of the resummed formula for the periodic orbit contribution to 
eigenfunctions obtained in the present work (5.34), and the analogous expression (2.24) Agam and 



Fishman |]13| obtained is very similar. Therefore, it is illuminating to compare the two, because the 



methods of derivation differ qualitatively. 

First, the expression (2.24) Agam and Fishman obtained in terms of the quantities that apply for 
tiling billiards with periodic boundary conditions on the pseudo-sphere will be presented. For this 
purpose, the quantities identified in Eq. 3.22a, 3.22b and the following, are required: 



1. g|| 2k as /c ^ oo; 

2. From Eq. 2.2C1| and the monodromy matrix for billiards with periodic boundary conditions (see 
Eq. 42 in H): i ' 



Rt-Rl 



k as k ^ oo; 



3. Up = Lp and 7^ = 0, because the simple case of periodic boundary conditions is being 
considered (see Eq. 3.20); 



4. From Eq. |zT| : C^^p{h,E) = kL 



_ Mi 
M.p 2 ■ 



Substituting these in (2.24), one obtains 

1 




ip Jkao'^, tanh ^i^L„ 



{p} M j=o ysinh(j + l)Lp 

Ak, 



ka 



2Q^{K,k^) 



L^^p + q± 



(7.1) 



3S kn 



00 



where {p} denotes the list of primitive periodic orbits of the billiard. We turn now to write the 
result (5.34), obtained in this work, in a form that can be compared with (7.1). The asymptotic 
expansion that led to (2.24) was obtained after it was argued that the leading order contribution 
was from a small region around periodic orbits |11]. Therefore, the relevant part of (5.34) in such 
a region should be found. It is easy to see that what remains of J2{p} J^her/Fp sum over the 
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isometries whose invariant geodesies constitute the actual trajectory in the billiard, because images 
of periodic orbits lie outside of the billiard, and therefore there are no points in their close vicinity. 
The remaining sum will be denoted by J2- The next stage is to obtain an approximate expression 



p.o. 



for the C-dependent part of (5.34). For djr{C) —>■ 0, it turns out that (from Eq. 3.6a) 



dgiC)^ Lg + tanhhgdf\c) 3S dj {() ^ . (7.2) 



This means that, to leading order: 



,jfca(dftj,j+i^-l(C)-a+l)ip) ifc<,tanhi(j+l)Lp ^.^{C) 



e 



Jsinhdhpj+ih-i (C) ^/sinh {j + l)Lp 

and 



(7.3a) 



In cosh 41(0 ~ d^liC) as d^,, (C) ^ . (7.3b) 
Inserting these into (5.34), one obtains: 



xe 



-i-KN{E)+ikaLfj,^p gj,£g 



2Q2(£a:„) (^'^'^ + '^^°-^^^ 2") 

as /cq ^ cx) and (i^o.(C) ^ • 



One can see that (7.1) and (7.4) are exactly the same, although the derivation of these two relations 
was completely different. 



The above comparison brings to light what is missing in the result of Agam and Fishman Their 
analysis was based on an expansion around periodic orbits, while in the derivation of (5.34) here, 
there was no such assumption. In particular, (5.34) includes the contribution of images of periodic 
orbits, which one would never expect to obtain from an expansion around periodic orbits, because 
these images tend to lie far away from points in the billiard. 

The expression (7.4) obtained near periodic orbits reveals the role of periodic orbits in the eigen- 
function density. As the distance from a periodic orbit is increased, fringes occur, that oscillate on 
a scale of and as /c ^ 00 they are washed out, as a result of contributions of various orbits. 
As in typical chaotic systems, one observes that the leading contribution should arise from points 
in the close vicinity of periodic orbits. From this analysis an analogy with typical chaotic systems 
can be drawn. As the structure of the formulae in those systems is the same, it can be expected 
that for some eigenfunctions only a small number of periodic orbits would contribute appreciably. 
Of-course this cannot be the common case, since usually many periodic orbits contribute, giving rise 



32 



to elaborate interference patterns. But, for the few eigenfunctions for which this is not the case, 
scars should be prominent. The question of the existence of scars in the systems studied here is 
addressed in the next section. 



8 A search for scars 



In this section, the possibility of finding scars, as usually understood by physicists, namely eigenfunc- 
tions that are appreciably large on unstable periodic orbits, will be explored The starting point are 
analytic calculations similar to those of for a billiard with periodic boundary conditions. Then, 
corresponding results are studied for a hard wall billiard. Finally these are tested numerically. 

Scarring of a wavefunction in a quantum system, means that the corresponding probability density 
is considerably larger than average near some unstable periodic orbit of the corresponding classical 



chaotic system. Agam and Fishman [15| suggested to quantify this concept and proposed to test 
whether an eigenfunction is scarred or not. They defined Yp{Ea) as the amount of wavefunction 
density in a narrow tube in phase space, with a large enough cross section around the orbit p, in 
excess of the mean value. The term "large enough" means that the tube has to be wide enough on 
the scale of oscillations as the distance from the periodic orbit is increased (governed by h). This 
tube cannot be too wide on the classical length scale if it is to encompass one short periodic orbit 
only. Using the semiclassical approximation for eigenfunctions (after the integral over the momentum 
coordinates is performed), the "scar weight" is given by (see AF 5.4): 

Y^{E^)= J df,iC)\i;^{C)\l , (8.1) 

tubep 

where |V'a(C)lpa is the contribution of the periodic orbits to the eigenfunction density. Eigenfunctions 
are normalized, therefore, if the number obtained is large (on a scale of 1) when compared to integrals 
over tubes around other periodic orbits, then a considerable fraction of the eigenfunction density lies 
in this tube and the eigenfunction should be scarred. If, on the other hand, the number is small, 
the eigenfunction should not be scarred. Agam and Fishman tested this criterion for a truncated 
hyperbola billiard. 

The effect of taking the narrow tube around the periodic orbit p is to remove from (5.34) all elements 
with large dp^ 's. There are two types of such elements: images of p, all other periodic orbits (and 
their images). This means that the tube removes from (5.34) all of the contributions that do 
not appear in (7.1). Therefore, the analytic results derived here concur with those of Agam and 
Fishman completely. Namely, if the tube is narrow enough, then the contributions from periodic 
orbits other than p would generally be negligible. The reason is that, because periodic orbits other 
than p typically have large dp^ , the phases in (5.34) (as can be seen clearly in Eq. 7.4) oscillate 
wildly. On the other hand, there is a narrow region around p in which dp^ is very small, so that the 
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phase in (5.34) is stationary, and therefore p is dominant. 

Instead of integrating over a complicated tube that encompasses many segments of geodesies in the 
billiard, the symmetry of the eigenfunction may be used to simplify it. Since the eigenfunction is 
invariant to group isometries it is the same in all the tiles obtained by them. Therefore, a tube of 
length Lp that extends from the billiard into nearby tiles may be constructed, just like one constructs 
a periodic orbit, in the following manner: 

1. Construct a tube around one of the periodic orbit segments in the billiard; 

2. Continue the construction into the next tile the segment continues into; 

3. Repeat this process in all the copies the orbit visits during one period. 



Therefore, restricting the integration to a narrow tube surrounding the invariant geodesic of p picks 
out of (5.34) the contribution of two transformations, p and p^^, provided that the tube is narrow 
enough. The reason for this is the behaviour of the function 

Ai?''''^(C, k) of Eq. 5.31c (here p' stands 
for any primitive element from the group) far away from invariant geodesic of hp'h^^ . This function 
involves a phase that causes it to oscillate wildly during the integration over the width of the tube. 
Only for the invariant geodesic of p (namely for the group elements p and p^^) can this phase be 
stationary throughout the length of the tube. This means that the contribution of p and p^^ is 
dominant. This argument is, of-course, the same argument for typical chaotic systems that Agam 
and Fishman |15| used. 



The resulting scar weight: 

Yp{E^)r.2 J dMO ^^J.^ A^^''HC,fej[ , (8.2) 



tuber, 



where p is the isometry whose invariant geodesic was chosen as the axis of the tube, and I is the 
identity. The factor 2 arises because the contribution of p"^ is exactly equal to that of p. 

It is always possible to use isometries to transform to coordinates in which p is diagonal. In such 
coordinates the invariant geodesic of p is the positive imaginary axis in the Poincare half-plane. 
Another isometric coordinate transformation may also be performed, so that the invariant geodesic 
section under consideration extends from 1 to e^p. 

Now (8.2) reads: 
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where: 



iy\/cosh«) — 1 

dy 



dx< 



ik[d^j+l{x+iy)~{j+l)Lp) 



-y\J cosh tu— 1 



X Erfc 



A; 



(8.4) 



where w is the width of the tube assumed narrow enough to capture only one periodic orbit. Note, 
that a tube around a geodesic is defined as the geometrical location of all points of constant perpen- 
dicular distance {d^) from it. In the coordinates in which p is diagonal, cosh^dp (C) is very simple 
(see Fig. D: 

cosh^ dp{C) = cosh^ ^CMCl ~ ^ + 
and depends only on the variable: 



-)■ 

yJ 



(8.5) 



yZ 



(8.6) 



Due to Eq. 3.6a, cosh d^^+i {() depends only on this variable too. (This result is correct for odd 
boosts too, due to Eq. |3.6b| .) 

As the tube is very narrow one may use the approximations ( [7.3a| ) and ( [7.3b| ). Changing variables 

X ^ X = - one obtains: 

y 



Lr. 



Vcoshui— 1 



sinh {j + l)Lp 
X Erfc 



J dxl e*'=tanh(2±lL,)x2 



2Q^K,k) 



\J cosh ui— 1 

L 



k^ v2 



as k 



oo 



(8.7) 



The integrand of (8.7) contains an exponent with a rapidly changing phase, so the method of 
stationary phase may be used. Since most of the contribution to such integrals comes from an 
exponentially narrow region around X = 0, then if the tube was wide enough on the oscillations' 
scale (■^), the integration limits may be taken as ibcxD without introducing a large error. In addition, 
the complementary error function hardly changes in the contributing region, so one may approximate 
it by it's value at X = 0, i.e. on the the periodic orbit. One then obtains: 



Lr 



sinh 



■Erfc 



k 



2Q^{K,k) 



as A: ^ cxD , 



(8.8) 
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where S,^,p{k) = (^L^,p — Inserting (|8!^ ) into (8.3), one obtains: 



3p 

xE 



Erfc 



as A;„ ^ oo 



(8.9) 



8.1 The weight of scars for hard wall billiards 



Using the the same definition for the scar weights ([SJ] ) that was used before, very similar results are 
obtained: 



2kaA'{Ea 



-i-KN{Ea)+ikaLfj,^ 



E 



osinh^L 



■Erfc 



P X 



C/i,p(^a) 



as ka —>■ (yo , (8.10) 



where 



self-retracing 
else 



(8.11) 



The only differences in the derivation and in the final result are numerical prefactors. For singular 
periodic orbits an additional factor of half is required, because they run along an edge and therefore 
only half of the tube around them is really in the billiard. This factor is compensated by the factor 
of two associated with the existence of an odd counterpart to every even singular primitive periodic 
orbit. The factors bp result from the fact that if periodic orbit is self-retracing, than p^^ is in the 
conjugacy class of p and therefore the prefactor two that appears in (8.2) is missing. 



8.2 A numerical search for scars 



In order to check for scars on some periodic orbits of the billiard defined in Fig. ||, the scar weight of 
Eq. 8.10 was calculated on the roots of the semiclassical approximation of the spectral determinant 
(see Section 33), that are the semiclassical approximation to the eigenenergies. The results for the 
four shortest primitive periodic orbits are presented in Fig. |7|. Also presented in the figure are the 
locations of the energies that admit standing waves on the respective periodic orbits. These are 
the energies satisfying the Bohr-Sommerfeld quantization rule. For billiards with Dirichlet boundary 
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conditions it reads: 



27r 



n+ -rih 



En 



(8.12) 



where Lp is the length of a periodic orbit and is the number of times it bounces off the walls 
during one period. 



The sum rule |15| 



1 



-TTp{Ea 



(8.13) 



was found to be satisfied to a very good approximation in every case studied. In Eq. 8.13 , p{Ea) 
is given by the derivative of N{Ea) (3.21) with respect to E. Eq. 8.12 is the expected result 
for ka —>■ oo, because the stationary phase limit was used to calculate the scar weight from it's 
definition, (|8.1D, i.e. integrating over a narrow tube and integrating over the whole billiard is the 
same as ka — > cxd. Using the fact that the eigenfunction is normalized, obtaining ( 8.13 ) is straight 
forward. 



The general features of the scar weight, Yp{Ea), on various periodic orbits for the eigenenergies as 
presented in Fig. |7| agree with the general features found in [|l^. Prominent among these is that the 
local maxima tend to be located at eigenenergies close to energies that admit standing waves, given 
by (|1]). 

The main difference between Fig. ^ and the figures presented in [|^, is that it seems that the 
functions found in the present work are not as smooth as in the previous work. As a result of this, it 
is difficult to find an eigenenergy for which Yp{Ea) is negligible for all p's but for a few, because the 
scar weights of various periodic orbit fluctuate strongly. The reason for the erraticness is not clear 
at this stage. As a result, it is difficult to implement Agam and Fishman's criterion for scars here. 

In order for an eigenfunction to be scarred it must have a pronounced contribution from only a few 
periodic orbits. Therefore the scarred eigenfunctions should be the ones for which the scar criterion is 



appreciable on one or two periodic orbits only. Agam and Fishman |15) found this to be in very good 
agreement with the features observed in the eigenfunction densities they calculated numerically, for 
a truncated hyperbola billiard. 

In addition to this, one must bear in mind that in order for a scar to be prominent it must be 
concentrated in a narrow tube around a periodic orbit, where "narrow" means narrow with respect 
to relevant classical linear dimensions in the billiard (for example, the distance between short periodic 
orbits). The width of the tube in which one may expect to see a scar is given by the scale on which 
the oscillations, as a function of the distance from the periodic orbit, occur. According to (7.4), this 
scale is approximately given by: 

ln.r ~ -i= . (8.14) 



kr) 
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The highest energy for which the results here were trustworthy (from the numerical calculation of 
A(£')) was approximately 

« 5000 , (8.15) 

that corresponds to 

lose ~ 0.12 . (8.16) 

For comparison, the distance between the point the shortest periodic orbit hits the side h and the 
point on that side where the second shortest periodic orbit hits, is 0.16 (see Fig. ^). This means 
that even at the highest eigenenergy for which the scar criterion could be applied with any reliability, 
within our numerical limitations, the tube where one expects to see scars is too wide for resolving 
individual periodic orbits. Nevertheless, a numerical search for scars was initiated. 



There are two ways to go about such a search. The first, which yielded good results in the past |15], 
is to try to use the scar weight graphs (see Fig. ^) to predict whether or not a prominent scar 
exists in some eigenfunction. The second method will be described below. The first was used here, 
without decisive results. When comparing the predictions with the figures of the numerically found 
eigenfunctions, sometimes scar-like features were observed where they were expected and sometimes 
they were not. Some figures of eigenfunction densities are presented in Fig. ^ and Fig. In principle 
these figures are slightly misleading, because the metric (3.1a) is not constant. But, as the billiard 
is not very extended, as can be seen in Fig. || and |8|, the aberration is not large, because: 



4r drd(j) 
(l-r2 



dl^= : =4(1 + 0[r^])r, 



As the shortest periodic orbits lie in a region with r <\ approximately, the deviation of the invariant 
measure from the Euclidean invariant measure is negligible. 

For example, according to Fig. |7| eigenfunction #84 should be scarred along the second shortest 
periodic orbit, while eigenfunction 7^68 and eigenfunction #90 should be scarred along the fourth 
shortest periodic orbit. In Fig. ^ one can see that indeed there are some high peaks of these 
eigenfunctions along the said periodic orbits. On the other hand, according to Fig. ^, eigenfunction 
#91 should not be very scarred along the fourth shortest periodic orbit., while in Fig. ^ one can see 
that it does have some high peaks along that periodic orbit. Another example is eigenfunction #92 
which, according to Fig. 0, should have exhibited a scar along the shortest periodic orbit and does 
not, as can be seen in Fig. |10[ In the same figure one can see that the next level, #93, does show 
a scar on that periodic orbit, as the scar criterion predicts. 



At this point it is worth mentioning that the scars observed in [|l5[|8| and more recently in [36| tended 
to be more pronounced than the scar-like features observed in the present work. 

One may also choose an opposite strategy in order to check the viability of the scar criterion. This 
strategy is faulty, as will be explained shortly. In this method one looks for scar-like features in 
eigenfunctions, and then checks the scar criterion graphs for a correspondence. There were several 
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discrepancies, one of which was mentioned above, namely eigenfunction #91 (see Fig. |^ which has 
a prominent scar-like feature on part of the fourth shortest periodic orbit. This strategy is misleading 
because there are many complicated structures in high energy eigenfunctions in a billiard, as they are 
the result of an elaborate diffraction pattern. This means that a scar-like structure may have nothing 
to do with a periodic orbit, as it may be just an accidental pattern. The first strategy, described 
above, filtered out some of these scar look-alikes automatically, because a large scar weight indicates 
that the structure observed may be related to a periodic orbit. 

The reason why one expects scars in the first place, despite claims to the contrary |17, 18,|l9|, is 



based on the observation that many features of the special class of systems under study here are 



similar to those of typical chaotic systems (for example Eq. 7.4). Since Agam and Fishman ||15| 
found scars in eigenfunctions whose energy admitted standing waves, one is led to expect the same 
here. As many prominent scars were not found in the system under numerical study, although the 
scar weight graphs suggest that scars should be present where standing waves occur, there may be 
important differences between the systems studied here and typical ones. 

An important difference that may be the cause for the scarcity of scars in the eigenfunctions examined 
is the very high topological entropy of billiards on the pseudo-sphere. For example, in the hyperbola 



billiard, a derivative of which Agam and Fishman |15| used, the number of periodic orbits up to 
length d is p7|| : 

O.Gd 

^(^)~777^ asd^oo. (8.17) 
0.6a 

This is much less than in tiling billiards on the pseudo-sphere, as can be seen in ( |5.22| ). This leads 
to the following consequences: 

1. It is difficult to find all the primitive periodic orbits up to some length; 

2. The periodic orbits are denser than in the typical case, so shorter wavelengths may be required 
to observe scars clearly. 



Evidently, it would have been advantageous to study states with higher energy, that corresponds 
to shorter wavelength, than was done in the present work. Unfortunately computational limitations 
prevent this, leaving the numerical results inconclusive, as to the existence of scars. 



9 Summary and conclusions 

In this work an exact expression for the eigenfunctions of special systems, whose classical dynamics are 
chaotic, was obtained. These systems, tiling billiards on the pseudo-sphere, are special in that many 
formulae, that for typical chaotic systems are correct only in the leading semiclassical approximation. 
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are exact for them. Specific differences between the exact expressions, as well as the approximate 
ones, and those used for general chaotic systems were identified explicitly. 



The main result of this work is Eq. 4.1(]| , an exact expression for (the diagonal part of) Green's 
function, in the form: 

whose residues are the squares of the exact eigenfunctions of the Hamiltonian ( |4.13| ). This result is a 
manifestation of Fredholm theory |Q for the systems studied here, because that theory predicts that 
Green's function in general can be expressed as such a ratio. The special feature of tiling billiards 
on the pseudo-sphere is that both numerator and denominator can be expressed in terms periodic 
orbits, without the introduction of any approximation. It would be interesting to study the exact 
connection between the results presented here and those of Fredholm theory. Such a connection 
would supply a recipe for reordering the problematic sums in the numerator and the denominator to 
ensure rapid convergence to the exact result. 



Summation over repetitions, leads to the form (|9.1| ) of the Green function and to the appearance of 
the spectral determinant in the denominator. Therefore repetitions directly determine the spectrum. 
Any specific zero is not related to a specific factor in the product of terms of the form of (|2.8|) , 
but rather to a zero of the the product that does not converge absolutely. This mechanism has 
been found in earlier work in the framework of the semiclassical approximation |13] and of Fredholm 



theory 10, but here it is established explicitly in the framework of the exact theory. 



It was shown that there are three different kinds of contributions to the exact J\f{C,E). The first of 
these is associated with the non-periodic orbits of the system, the second with the periodic orbits, 
and the third with their images. The first two contributions appear in the leading semiclassical 
approximation for Green's function (and for eigenfunctions) of typical chaotic systems. 

The third contribution, that of images of periodic orbits, can be given the following interpretation. 
Such a contribution can be associated with a non-periodic closed orbit, whose length is c^/iph"i(C)- 
In typical chaotic systems this contribution is approximated by nearby periodic orbits, but in a tiling 
billiard on the pseudo-sphere, the fact that hph^^ belongs to a conjugacy class that corresponds to 
a specific periodic orbit, provides an exact recipe for accounting for it's contribution. In general, such 
a correspondence cannot be found with the standard semiclassical approximation methods, because 
the closed orbit hph^^ is associated, with need not lie close to the periodic orbit to which p is 
related. 



In the general result obtained for chaotic systems, contributions are taken only from the close vicinity 
of periodic orbits, assuming that contributions from distant ones cancel due to destructive interference 
|jlO|,|ll|. Here, such an assumption was not required, since the contributions from all periodic orbits 
can be calculated exactly. If further approximations, namely an expansion around periodic orbits, are 
made, the expressions reduce to those found for general chaotic systems. For general chaotic systems. 
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the results can be obtained only in the framework of such an approximation. For the systems that 
are studied here, the only correction to the semiclassical results are the contributions from images of 
periodic orbits, namely orbits outside the billiard. 

The second result obtained is a semiclassically approximate expression, resulting from an effective 
truncation of the infinite sums over periodic orbits of the exact J\f{C,E) and A(E') by the Berry- 
Keating resummation 0. The truncated expressions were compared with previous results for typical 
chaotic systems IQ. As the previous results were derived assuming that the dominant contributions 
arise from the close vicinity of periodic orbits, the same limit was checked for the results obtained 
here. The two expressions turned out to be identical. 

The resummed expression for eigenfunctions of hard wall billiards was checked numerically for the 
two lowest states, and reasonable, although not perfect, results were obtained. In principle, higher 
orders of the semiclassical approximation can be derived for the systems studied. Inclusion of these 
in the numerics should yield better results for the low eigenstates checked. 

In addition, a criterion for scars ||l^ was derived for the chaotic systems at hand. It was tested 
numerically to check whether scars exist in the systems studied. Inconclusive results were obtained, 
as the energy attainable in this work was not high enough, due to computational limitations. In order 
to resolve the question of scars, better numerics are required. One way to go about this would be 
to work with arithmetic billiards, although they are not generic, as for them there exist systematic 
methods for obtaining long lists of periodic orbit lengths, that are complete. This should enable the 
study of states with high energies. 

Throughout this work the nature of the exact results obtained for tiling billiards on the pseudo- 
sphere was found to be very similar to results that are only semiclassical approximations for typical 
chaotic systems. Therefore, tiling billiards on the pseudo-sphere may be useful as a testing ground 
for approximate results concerning typical chaotic systems. 
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Appendix A The derivation of Eq. 2.22 



In this appendix the sum over m in en.^^"'^ Yl, cjf'"^^g„ 



of (2.15) will be performed, 



leading to (2.22). Since the sum over fi and m converges absolutely, the order of terms can be 
changed. In particular, exchanging the summation over composite orbits with the summation over 
m and using Eq. 2.19b, one finds: 



m=0 



— XRf^X 



oo 

X ^ 9m 

m=0 



— XR^X 



jp 

j=0 



Exchanging summation order again, one obtains: 



m=0 



— XR^X 



3p 



j=0 



m=0 



— XRf-,X 



giXRpXg-m(j + l)Mg: 



Comparing the curly bracketed term with (2.21) it is easy to see that 

exp 



m=0 



— XRf^X 



en 



|XRpX -m(j + l)Uj: 



■ ~ AfJ + ^ T 



so that in conclusion, since in two dimensions (see AF 3.11 ii 3.13): 



exp 



exp 



^XRpX tanh ^-^u^^ 



det(M^+^ + /) 



2 cosh ^^tif 



one obtains 
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m=0 



-XRpX — ^^~p\ 

j=0 



exp 



I^XRpX tanh ^-^Up 



2 cosh ^-^Ug^ 



E 



i=o 2 cosh iii-up 



exp 



— XRqX tanh ^^-^^^ — n„ 



C^^'^^ was defined in (|Z2^ ). Substitution of the last equation in (2.15) leads to (2.22). 



Appendix B Writing Legendre functions as sums over e 



To prove (4.2) one starts with the integral representation of the relevant Legendre function [38| 



Q-l-ik [coshd] = -7^ J 



y/2 J ^cosh(d + y) — cosh((i) 



for d > 



which can be rewritten: 



n=0 



The summation and integration order can be reversed, because the integral and the sum are absolutely 
convergent: 



Q.,, [coshd] = E ^^^;^e^'-^-'"^' I dy 



,-{n-ik)y 



The last integral is known and one obtains 



Q-HJcoshd] = Ei^y^i? 



1 1 



g{ifc-i-2n)(i 



(B.l) 
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Since the d's to be used in the above expression are the dg{C)'5 of Eq. p.6a| , one can write, after 
some straight-forward algebra: 



where 



N2l 



1 + W1 + 



+ 



cosh^ d^ (C) (1 - a;)2 cosh^ d^ (C) (1 - x) 



Now 



ik-i-2n 



can be expanded in powers of x = e 



cosh2dJ-(C) 



/ x,dUo 



m=0 



1 , 1 

Ak 2n 



cosh^d^iC) 2 



where the first few am[x,y]'s are: 



oo [x, y] = 1 , 



(B.2) 



Since the ajn[x,y\s are polynomials of y, one may write: 

am[x,y\ = am[x,dr^]e^y 



77=0 



to bring (B.2) to a more convenient form: 

1 



m=0 



cosh2d^(C) 



cosh2d^(C)^ 



ikLg -(\+m+2n)Lg 



X e ''e ^2 

??=0 



. (B.3) 



Inserting (B.3) into (B.l) one obtains: 



OO OO /r) ^\\\ 

Q_i [coshd,(C)] = E E ^^^^^ 



n=0 m=0 
1 



cosh2 4(C) 



(2n)!! 
drj 



1 1 
n H zfc, — 

2 ' 2 



X 

i+2?i-ifc 



^cosh2(iJ-(C)^ 



*^-^9p-(|+'"+2ri)L9 



xe "e ^2 

17 = 



(B.4) 
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Changing summation variables in (B.4): 



\ ^ \ ^ . „ m'=m+2n 
n=0 m=0 

leads to Eq. 4.2. 

The alterations for odd boosts 



oo oo 



m' ~2n 

71=0 ■m'=2n 



E E ^nBm'-2n , 



(B.5) 



m'=0 n=o 



In deriving ( |4l2|) , (|l6i ) was relied upon to obtain an expansion of e'^"^''^ in powers of e ^p, when 
the group included only even boosts. For odd boosts ( p.6b| ) will be used, leading to: 



g^odd 



x(l + x) 



1 + Jl 



X 



cosh^ 4 (C) (1 + x)2 cosh^ 4 (0(1 + 



(B.6) 



One can see that the only difference between (B.6) and (B.2) is the change x — > —x. This means 
that when (B.6) is expanded in powers of x, an expansion similar to (B.2) will be obtained, but for 
a factor (—1)™. Therefore: 



g odd 



The last equation means that the only alteration that (4.2) requires for odd boosts is the multipli- 
cation of each term in the series by (—1)'". Eq. 4.15 and 4.16 follow from this immediately. 



Appendix C The coefficients C„ and Cjf'"^^ 



The coefficients are defined in ( 2.10b| ). For tiling billiards on the pseudo-sphere they are given 
by: 



-2L. 



(C.l) 



where the definitions in (3.20) were used, and the function dj^ {x) was defined in ( |2.19c ). Notice 



that the C^'s for the billiards under consideration are real numbers. 
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The coefficients cjf'"^^ (p stands for p, q or s) are also all real. Using their definition of (|2.19b ) 
and the definitions of Eq. 3.20, they are given by: 



j=0 

3q 

i=o 

is 

= C^_,(-lpe-t(^'^+i)^^ ^(-iye-2™(j+^)^»dj,_j . (3.4c) 



Appendix D Resummation of the numerator of Eq. |4.iU| 



It is useful to write Npo{C,E) of ( |4.12D in a form where the dependence on k is shown explicitly: 

1 



-1 ^ L J 

{p} m=o/ier/rp n=o 



27r 



1 1 

n H ik,— 

2 ' 2 



(D.l) 



rj=0 



where 



and the definition of jr„ 



i;i,p(C, ^) = + In (cosh" f^hph-i (0 j + 
A;,d^_ip^(C) of (4.3) was used to define: 



(2n-l)!! 



h2ra 



2/. ^] = E (2nl!! ^] 

n=0 ^ 

Using ( pTTl] ) and (D.l) in (5.28) one obtains: 

Ar(£;(fe))^-i^^[[/^(fc) + C/^(-fc)] , 



(D.2) 



(D.3a) 



1 oo 

A/-p„(C,£;(fe)) 5] EE[^M,™,M(C,fc) + C/M,m,/.(C,-A:)] , (D.3b) 



{p} /iGr/rp»n=o ^l 
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where: 



-oo+i(l/2+e) 



oo+i(l/2+e) 

(1) e-^^{^iEiz+k))-N{Eik)))+izLf, 



n=0 



,dri,ri 



-i7rJV(£;{fc))+ifcL^,p(C,»?) 



-oo+i(l/2+£) 

^^-iiT{N{E[z+k))-N{E{k)))+izLl,^j,{C,ri) 



oo+j(l/2+£) 



1 ./ ,^ 1' 
n + --^[z + k),- 



\t turns out that: 



u,{-k) = u;{k) , 



This is because the C^'s and cjF'^^^'s are all real (see [Appendix"^ ). In addition ||7| 



N{E{-k)) = -N{Eik)) , 



and !§: 



oc 



dt e 



ikt 



y/2J Vcosht - 1 





B[n + 1/2 -ik, 1/2] 



jji~l/2—ik 

dt = B*[n + 1/2 + ik, 1/2] . 

vl — t 



Inserting (D.5a) and (D.5b) into (D.3a) and (D.3b) yields the manifestly real expressions 

M{E{k))=-^Y.'^{u,{k)} , 
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Mpo{C, E{k)) = — E E EE^ ^p,w(C, k) . (D.8b) 
^ {p} her/Fp m=o fj. ^ ^ 

The next aim is to obtain an approximation for the integrals in (D.4a) and (D.4b). From Weyl's 
asymptotic series ( |3.16| ) one obtains: 



vr {N{E{z + k)) - N{E{k))) ~ —z + ^z'^ as k ^ oo 



Ak 

Y 



(D.9) 



Notice that the series truncates after because the next terms are all exponentially small in k |2C]. 
This fact is specific to billiards with periodic boundary conditions, and will not hold for billiards with 
reflecting walls. 

Using the integral representations of the special functions in the integrands, (D.7a) and (D.7b), 
double integrals are encountered. Since the integrals are all absolutely convergent, the order of 
integration can changed in each case. Then one encounters, in the leading order, an integral that 
can be evaluated analytically: 



lim 



-oo+i(l/2+e) 



dz _qi. 
— e 2k ' 
z 



1 



Erfc 



oo+i(l/2+e) 



'A 

2Q2 



(D.IO) 



Using ( |D.10D in Eq. D.4a and D.4b one finds: 



dy e^^y 1 
\/2 J V cosh y — 1 2 



Erfc 



k 



2Q^K,k) 



,(D.lla) 



n=0 



,d7],ri 



n-l/2-ik -j^ 

X / dy^^—== — Erfc 



2Q2 



(^(C(^,C,^)-ln(y) 



as A; ^ oo , 



77 = 



where: 



h A 



(D.12) 
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To leading order, (D.lla) is given by: 



2 ^ 



k 



2QHK,k) 



, as A; ^ cxD , (D.13) 



which yields Eq. |5.30b when inserted into Eq. D.8a. Similarly, the leading oder of (D.llb) is given 
by: 



1 



n=0 



, drj, rj 



xB 



n -\ ik, — 

2 2 



■Erfc 



k 



2Q2(K^ k)^i^,P 



r)=0 

as ^ oo . (D.14) 



In (D.14) the 9^ differentiations must now be performed. But since the objective here is only 
to derive a leading order approximation, the derivatives of the complementary error function can be 
neglected. The reason for this is that the derivatives of complementary error functions are proportional 
to Gaussians centered their cut-off. Since their contribution is proportional to their width. The 
Gaussians' contribution is smaller than the leading term, proportional to the complementary error 



function. Thus, the leading order approximation of ( 4.12 ) is Eq. 5.31a 



M,o{C,Eik))^^^J2 E ^^'^'\C,k)\ ask 
^ I {p} her/Fj, J 



oo 



where, using (D.2): 



A'-'(C..).-.£gi^B 



1 ., 1 

n-\ IK,— 

2 ' 2 



1 J, 1 o 



2QHK,k) 



(D.15) 



In the last formula, the definition: C/j,p(^)C) = C|!,p(^)Ci^ = 0) was used. 

Using arguments similar to those leading from Eq. 2.15 to (2.22), Eq. D.15 can be simplified to yield 
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Eq. 5.31b (for the details see |Appencjix~^ ): 
where: 



coshd;,pj+i^_i(C) 



(D.16) 



The Ci^'^')' s were defined in ( [2.23| ) and the definitions of (3.20) for even boosts were used. 



In order to be consistent the leading asymptotic behavior of the functions appearing in Eq. D.16 
must be used. Since (BV G.21): 



one obtains Eq. 5.31c. 



Appendix E The derivation of Eq. D.16 

In this appendix, Eq. D.16 is derived from Eq. D.15. Since A^P''^\() involves only absolutely con- 
vergent sums, the summation order can be exchanged and the sum over /i performed last. Using the 
definition of ^fj (A:,C) (after Eq. 5.31b) and the definition of A|f (D.2) in Eq. D.15, one obtains: 



A(P''^)(C, k) = Y, ^'f''"HC, k)e-'''^^^^+'^^^''^ Erfc 



2Q^iK,k) 



(E.l) 



where 



m=0 



n=0 



(2n)!! 



1 ., 1 

n H ik,— 

2 '2 



1 1 
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Using the definition of the C^^ '^^'s (see Eq. |2.19b| ) yields 



3p 



j=0 



where: 



m=0 



n=0 
1 



1 ., 1 

n H IK,— 

2 '2 



,ik - 1/2 - 2n 



^+2n-ik 



Exchanging the sum over m with the sum over n and shifting summation (opposite to the operation in 
Eq. B.5) yields an expression similar to the one that appears in (B.4), that led us to (4.2). Therefore: 



Tj [Lp,k] = e 



-ifc(j+i)Lpgi(j+i)ip 



ik 



cosh d 



(C) 



Insertion of the last expression into Eq. E.2 and E.l with the definitions of Eq. 2.23 and 3.20, leads 
to Eq. D.16. 
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Figure captions: 



Figure 1: Left - The action of some isometries of the group generated by reflections at the sides 
of the shaded triangle on it. One can see that the copies tile the pseudo-sphere. Right - The 
billiard chosen for the numerical calculations (a magnification of the shadowed area on the left). 



Figure 2: d^(C) is the distance between a point C and the invariant geodesic of g measured along 
a geodesic that is perpendicular to it. The invariant geodesic of g is aligned along the y axis. 



Figure 3: Comparison between the zeros of A(£^) approximated by ( 2.11| ) and the exact eigenen- 
ergies (empty circles) found numerically. 



Figure 4: The integration contours used for resummation. 

Figure 5: Comparison of an exact cut through the ground-state {Eq ~ 118.9) density (thin line) 
with the semiclassical result, before including periodic orbits (dashed line) and after including 
them (thick line). On the right is the exact ground-state density in the whole billiard, with the 
cut through it denoted by a line. 

Figure 6: Same as Fig. |5|but for the first excited state {Ei « 201.1). 

Figure 7: The scar weight on the four shortest primitive periodic orbits {p = 1 denotes the short- 
est, p = 2 the next shortest, etc.) calculated on the zeros of the approximate spectral determi- 
nant. The vertical lines are the locations of the eigenenergies predicted by the Bohr-Sommerfeld 
quantization rule ( ^.12| ). 

Figure 8: The four shortest periodic orbits: solid line - the shortest orbit, dashed - next in 
length, dot-dashed - longer still and the fourth shortest is dotted. 



Figure 9: (a) - Density of eigenfunction #84 {k = 67.0) and the second shortest periodic orbit. 
The fourth shortest periodic orbit on the density of eigenfunction ^68 (k = 60.8) (b), #90 
{k = 69.1) (c) and #91 {k = 69.5) (d). 



Figure 10: Density of eigenfunction #92 {k = 69.9) (a), and density of eigenfunction #93 
{k = 70.2) (b). On both figures the shortest periodic orbit is presented, as well as a tube of 
width around it. 
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